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Abstract 

The Hylleraas coordinates s = t\ + r2, t = r\ — r2, u = |ri — r2 1 are the natural coordinates for 
the determination of properties of the Helium atom, the positive ions of its isoelectronic sequence, 
and the negative Hydrogen ion. In this paper, we derive a new expression for integrals representing 
properties such as the energy, normalization and expectation of arbitrary operators, as written 
in the (s,t,u) coordinates. The expression derived is valid for both finite and infinite space. 
The integrals for the various properties are comprised in each case of two components A and 
B. The contribution of these components to the volume of integration and the normalization 
of a wave function for finite space, and in variational calculations of the ground state energy of 
the Helium atom confined in a finite volume is demonstrated by example. We prove that when 
the integration space is infinite, the expression for the energy and other properties employed by 
Hylleraas corresponds only to that of integral A. We further prove that for the approximate 
variational wave functions employed by Hylleraas and other authors, the contribution of the term 
B vanishes. This contribution also vanishes for the exact wave function. It is interesting to note 
that the component B to the integral is not mentioned in the literature. A principle purpose of 
the paper, therefore, is to point out the existence of this term. 

PACS numbers: 
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I. INTRODUCTION 

The Schrodinger equation for the Helium atom, or the positive ions of its isoelectronic 
sequence, or of the negative Hydrogen ion, is 

HV = EV, (1) 

where the Hamiltonian H is 

+ (2) 
2 2 n r 2 r 12 

where ri, r 2 are the coordinates of the two electrons, ri 2 is the distance between them(see 
Fig.l), Z the charge of the nucleus, and \1/ and E the eigenfunction and eigenenergy of 
the system. As noted by Hylleraas [2j, the Schrodinger equation in this instance is not a 
six-dimensional but rather a three-dimensional eigenvalue problem. The natural coordinates 
of the system are r 1; r 2 , and r 12 forming the triangle (r 1; r 2 ,r 12 ). The other three angle 
variables which determine the orientation of the triangle in space, because of the symme- 
try of the system, are arbitrary and can be integrated over. Equivalently, the system is 
uniquely described by the coordinates r 1; r 2 , and # 12 , the angle between the vectors and r 2 . 

As is well known, there is no closed-form analytical solution to the two-electron eigenvalue 
problem. Hence, the energy is most accurately obtained by variational minimization of the 
energy functional E[^f] defined as 

< 3 > 

with respect to the parameters in the wave function. In performing such variational calcu- 
lations, Hylleraas further introduced the 'elliptical' coordinates 

s — n +r 2) t — ri - r 2 , u— ri 2 . (4) 

Then assuming approximate wave functions dependent on these coordinates, he applied the 



variational principle to obtain the energy. T 



nere 



(s, t, u) coordinates over the past decades 



mm 



ias been similar work employing these 



In this paper, we investigate integrals of the form J fdr appearing in Eq. (3) in the 
elliptical coordinates (s, t, u) for both finite and infinite space. We have discovered another 
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way of writing the integrals in these coordinates that differs from the commonly employed 
one in the literature. In section II we derive the integral to be a sum of two contributions 
A and B valid for both finite and infinite space. This is a natural separation of the 
integration domain. We demonstrate the correctness of our derivation by three examples 
for finite space as described in section III. These examples correspond to (a) volume of 
integration, (b) normalization of a wave function, and (c) variational calculations of the 
ground state energy of the Helium atom confined to a finite volume. In section IV, we 
prove for the case when the space is infinite that the energy and normalization integrals 
employed by Hylleraas and others corresponds only to that of integral A. We further prove 
that for the approximate wave functions employed by these authors, the contribution of B 
vanishes. As the asymptotic structure of the exact wave function is known to also decay 
exponentially [9], once again the contribution of integral B for infinite space vanishes for the 
true wave function. We summarize our results and conclusions in section V. 



II. INTEGRALS IN (s, t, u) COORDINATES 

Let us first focus on the volume element dr of Eq.(3). Employing the symmetry of the 
system, dr can be rewritten as 

dr = dr 1 dr 2 = J(r 1 , r 2 , 8 12 )dr 1 dr 2 d6 12 , (5) 

where J(r±, r 2 , 8i 2 ) is the Jacobian of the coordinate transformation. By fixing ri first, carry 
out the integral of r 2 with respect to r x to arrive at (see Fig.l) 

J*r = J dr^r 2 = 2. / A, (6) 

Next, by integrating over r 1; we note there is no dependence on the Euler angles, so that 
Eq.(6) reduces to 

J dr = 8n 2 J r\r\dr\dr 2 sin9i 2 d9 12 . (7) 



We now introduce the new variable u, the distance between ri and r 2 : 

u 2 = r\ 2 = r\ + r\ — 2r\r 2 cos9\ 2 . 
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(8) 



FIG. 1: The coordinate system indicating the positions rx and r 2 of the two electrons with the 
nucleus at the origin. 



During the previous integration steps, if r\ and r 2 are fixed first, then 

u du = r x r 2 sin9 12 d9 12 , (9) 
so that (for the case of infinite volume) 

poo poo pir poo poo /*7*i+r2 

dr = 8n 2 / / r\r\drxdr 2 / sind^dOn = 8ir 2 / / rir 2 dridr 2 / udu. 

JO JO Jo Jo Jo J\ri-r 2 \ 

Eq.(lO) appears in the work of Hylleraas[2j. 

Let us next confine the two electrons in some finite volume of space such that < r\ < R 
and < r 2 < R. The reason we choose the same R is because of the symmetry between the 
two electrons. The limit R — > oo then leads to the infinite space integral. For the elliptical 
coordinates s and t of Eq.(4), it is easy to show that 

dr\dr 2 = -dsdt, r\r 2 = -{s 2 — t 2 ). (11) 

_ 4r 



FIG. 2: The domain S of integration of the coordinates 7*1 and r 2 . 



Therefore, for any integrand /, which is a function of (s,t,u), we have from Eqs.(lO) and 
(11) that the integral in terms of the new coordinates (s,t,u) is 

pR rR rri+r2 C C 1 f s S 2 — t 2 

fdT = 8ir 2 / rir 2 dr 1 dr 2 / uf(s,t,u)du = 8n 2 / -dsdt / f(s,t,u) — udu. 

Jo Jo J\n-r 2 \ J J 5' 2 J\t\ 4 

(12) 

Here S' denotes the integration domain in the new coordinates (s,t). This do- 
main has changed from S : [0 < T\ < R, < r 2 < -R] as shown in Fig. 2 to 
5': [0<s<2 J R;- J R<t< J R]. 

The function f(s, t, u) could be ipHip or if) 2 . As stated by BetheQ): " The exact symme- 
try requirement [of the wave function] then takes the simple form that ip be an even function 
of t for parahelium, an odd function of t for ortho-helium. Since the Hamiltonian is an even 
function of t and since the integrals in Eq.(3) contain two factors, the contribution to the 
integral from —t is identical with that from +t. We therefore restrict ourselves to positive 
values of t in the integrals and multiply the volume element by a factor of 2" . With that in 
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FIG. 3: The domain S' of integration of the Hylleraas coordinates s and t. 



mind, Eq. (12) can be rewritten as 

/ fdr = 2tt 2 [ [ dsdt [ f(s,t,u)(s 2 -t 2 )udu 
J J J s> J\t\ 



S> J\t\ 
-R ps 



= 2tt 2 / ds dt f(s,t,u)(s 2 -t 2 )udu 
Jo Jo Jt 

p2R p2R-s ps 

+ 2tt 2 j ds j dt j f(s,t,u)(s 2 -t 2 )udu. (13) 
Jr Jo Jt 



Now we can see from equation (13) that the integral has two components. For conve- 
nience, let us denote them as A and B where 

pR ps ps 

A = 2n 2 ds dt f(s,t,u)(s 2 -t 2 )udu, (14) 
Jo Jo Jt 

and 

p2R n2R-s ps 

B = 2n 2 ds dt /(s,t,w)(s 2 - t 2 )udu. (15) 

Jr Jo Jt 



We next give examples to demonstrate the correctness of this separation of the volume 
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of integration. 



III. EXAMPLES 

(a) Volume of integration 

We first consider the volume of integration by choosing f(s,t,u) = 1. The result should 
be the square of the volume of a sphere with radius R. Indeed, one easier way to determine 
this is by employing Eq.(7) 

jdr = 8n 2 -2(^-) 2 = (^f) 2 . (16) 

On the other hand, we can carry out the integral via its components derived in Eq.(13). On 
using the following integrals 

rs (2_ f 2\2 

J (s 2 -t 2 )udu= [ -^^, (17) 
d^—^- = - 45". (18) 



[ ds [ 

Jo Jo 



2R r 2R-s /„ 2 +2\2 00 d6 



f f 

/ ^ 

J R JO 



, (s 2 -t 2 ) 2 38R b 

dU —Y^ = ~45~' (19) 



we also have 



dr = A + B = 2tt 2 ( — + -^-) = (— ) 2 - (20) 

The agreement of the two ways to obtain the volume demonstrates that the new way of 
expressing the integrals within the (s,t,u) coordinates as defined by Eq.(13) is correct. 



(b) Normalization 

Let us next consider the normalization of the trial wave function 

if, = Ce~ as , (21) 

where a is a variational parameter, C is the normalization constant. Then in spherical polar 
coordinates for a finite volume of radius R, the normalization integral is 

1 = J i?dr = C 2 (4tt j\~ 2ar r 2 dr) 2 

= ^[l-e- 2aR (l + 2aR + 2a 2 R 2 )} 2 . (22) 



For R — > oo, C = a 3 /it, the well known result. 

In the (s,t,u) coordinates, the contribution from A of Eq.(14) is 

A = f i; 2 dr = 2tt 2 C 2 [ dse~ 2as f dt f duu{s 2 - t 2 ) (23) 
J a Jo Jo Jt 

(J2 2 -2aR 

= — — [1 (15 + 30aR + 30a 2 R 2 + 20a 3 R 3 + 10a 4 i? 4 + 4a 5 R 5 )}. (24) 

cr 15 

The contribution to normalization from B of Eq.(15) is 

r- f2R p2R-s ps 

B = 4> 2 dr = 2ix 2 C 2 / dse~ 2as / dt duu(s 2 - 1 2 ) (25) 
Jb Jr Jo Jt 

= C L^[ e -^R{i + 2 aR + 2a 2 R 2 ) 2 
or 

e -2aR 

rz-(15 + 30ai? + 30a 2 R 2 - 20a 3 R 3 - 10a 4 R 4 - Aa 5 R 5 )}. (26) 

It is easily verified that the sum of A and B of Eq.(24) and Eq.(26) are equivalent to that 
of Eq.(22). This again demonstrates the correctness of the two integrals derived. 

(c) Variational Calculations 

In the final example, we consider the Helium atom confined within shells of volume of 
integration corresponding to R = 2.0,3.0,4.0, and 5.0 (a.u.). We determine variationally 
the ground state energy of the atom thus confined. The approximate wave function we 
employ is that of Eq.(21), which then does not satisfy the boundary condition of vanishing 
at the surface. The total energy of Eq.(3) is determined in both spherical polar coordinates 
as well as via Eq.(14) and Eq.(15) of contributions A and B. A point to note is that the 
commonly employed expression for the integrand of the kinetic energy in (s,t,u) coordinates 
is valid only for infinite space and assumes that the wave function vanishes there. However, 
since the volume of integration is finite, there is a contribution to the kinetic energy from 
the surface of the volume as the wave function does not vanish on the surface (see the 
Appendix for details). The analytical expression for the kinetic energy contributions from 
integrals A and B and the surface contribution as well as those for the total energy are given 
in the Appendix. The various numerical contributions of the kinetic energy T, external E e 
and electron-interaction E ee potential energies are listed in Table I. The well known[] 
variational result for infinite space quoted to six significant figures is E = —2.84765 (a.u.) 



Ep^t 
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TABLE I: The variationally determined ground state energy of the Helium atom and of the corre- 
sponding kinetic energy T, external E ex t and electron-interaction E ee potential energy components 
as a function of the radius R of the volume of integration in (a.u.). The contribution from the 
integrals A and B for each energy component as well as the surface contribution to the kinetic 
energy are quoted as are the energy minimized values of the variational parameter a. 



R 


a 


T 


E ex t 




E 






A 


B 


S 


A 


B 


A 


B 




2.0 


1.652 


1.91327 


0.81583 


0.16747 


-5.64802 


-1.16273 


-0.88250 


0.20120 


-2.83048 


3.0 


1.687 


2.68121 


0.16476 


0.01175 


-6.59864 


-0.16328 


1.03104 


0.02713 


-2.84603 


4.0 


1.688 


2.82820 


0.02114 


0.00071 


-6.73641 


-0.01643 


1.05256 


0.00267 


-2.84756 


5.0 


1.688 


2.84729 


0.00205 


0.00004 


-6.75073 


-0.00131 


1.05480 


0.00021 


-2.84765 



for a = 27/16. For each value of R, the sum of the contributions from integrals A and B for 
the separate kinetic plus surface S contribution, external, and electron-interaction energy 
components is the same as those obtained in spherical polar coordinates. 

Note there is a contribution to each component of the energy from both integral A and 
B for each value of R. Due to the exponential decay of the wave function, the contribution 
from B diminishes but always is finite even at R = 5.0 (a.u.). Observe that at R = 5.0 
(a.u.), the total energy is the same as the well known result to six significant figures. 
However, the contribution of B becomes negligible in the R — > oo limit. Thus, it is only in 
this limit that the integral corresponding to B can be ignored. The variational results once 
again prove the correctness of the derivations for the two integrals A and B. 



IV. COMPARISON WITH THE HYLLERAAS INTEGRAL 

In this section, we compare our new way of writing the integral with that of Hylleraas. 
As noted previously, Eq. (10) appears in the work of Hylleraas. Following this equation, 
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and without any further detail, he then writes the integral for the energy and normalization 
for when the space is infinite as 

fdr = 2n 2 / ds du dt u (s 2 - t 2 ) f(s,t,u). (27) 
Jo Jo Jo 

Observe that in this integral, the integration over the variable t is performed before that 
of the variable u. (We note that in all subsequent literature employing these elliptical 
coordinates, it is the Hylleraas expression that is employed.) This order of integration is 
surprising because the variable u depends on the variables s and t. Hence, it is natural to 
perform the integral over u prior to that of t and s as in our expression Eq. (13). The 
derivation of Eq. (13) is a consequence of our attempt to understand how Hylleraas arrived 
at his expression of Eq. (27). For a comparison of Eqs. (13) and (27), we take the R — > oo 
limit of Eq.(13). The two results ought to be equivalent. However, we find that it is only 
UvriR yoo A (see Eq.(14)) that is equivalent to Eq.(27) of Hylleraas. For R — > oo 

POO J*S PS 

A\ R ^ 00 = 2ir 2 / ds dt duu (s 2 -t 2 ) f(s,t,u). (28) 
Jo Jo Jt 

As shown by Fig.4, in Eq.(27), for each u with < u < s, we have < t < u , the 
area swept is the upper shaded triangle. On the other hand, in Eq.(28), for each t with 
< t < s, we have t < u < s , so that the area is the same as in Eq.(27). Therefore, for 
any integrand, the integrals of Eq.(27) and (28) are the same. Thus, the Hylleraas integral 
is the same as that of A for R — > oo. The difference in the order of integration between 
the expression of Hylleraas and of our derivation is of critical significance for work on the 
Helium atom to be published elsewhere. 

Let us now turn to Eq.(15) and examine the integral B in the same limit. As we have 
shown earlier, in finite space, the contribution of B plays an important role. When R — > oo, 
we have 

BIr^ = 2vr 2 /im i? __ >00 / ds dt f(s,t,u)(s 2 -t 2 )udu (29) 

Jr Jo Jt 

r2R 

= UmR^^ 2n 2 / h(R, s)ds (30) 
Jr 

It is only when the function h(R, s) decays in a manner such that lim R >QO B vanishes that 

the Hylleraas expression corresponds to the exact value of the integral J fdr. Thus, for 
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FIG. 4: The domain of integration of the Hylleraas coordinates t and u. 



example, if one assumes as in all prior literature, that h(R, s) decays exponentially, and is 
of the form 

h(R } s)=g(R)e- as Y,*\ (31) 

l>0 

then since for any non-negative polynomial of s, 

limR^oo / dse~ as s l = 0, (32) 
Jr 

we have B\r >OQ = 0. However, if the dependence on the coordinate s contains terms of 

the form s p with p > —1 or is multiplied by some power of Ins, the contribution of B 
could be nonzero. For example, the integral J R h(s)ds for h(s) = 1/s is ln2 irrespective 
of how large R is. As noted above, for the form of approximate wave functions employed 
in the literature, the contribution of B always vanishes in the limit as R — > oo. It also 
vanishes for the exact wave function whose asymptotic decay is known [l| to be r /3 e _Qr , 
where (1 + 0) = (Z — N + l)/a, and a = \/2I, with / being the ionization potential. 

It is interesting that the contribution of the integral B to the integral J fdr is not 
mentioned at all in the literature. One would imagine that in spite of the fact that the 
contribution of B vanishes for the choice of a particular approximate wave function, mention 
of this term pointing out its lack of contribution would appear somewhere. But this is 
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not the case. Hence, one principal purpose of this paper is to note the existence of this term. 



V. CONCLUSION 

In this paper we have derived a new expression for integrals representing the energy, 
normalization, or other expectation values when written in the Hylleraas (s, t, u) coordi- 
nates, i.e. J f(s, t, u)dr. The expression derived for the intergral is valid for both finite and 
infinite space. The integrals for the various properties are comprised in each case of two 
components A and B. The contribution of these components to the volume of integration, 
the normalization of a wave function, and in variational calculations of the energy, for finite 
space is demonstrated by example. We prove that when the space of integration is infinite, 
the expression used by Hylleraas and others corresponds only to that of integral A. In the 
literature, the form of the approximate wave functions are such that the integrands for the 
energy and normalization are usually of the form e~ as Y^i m n hi(s)t m u n ; hi(s) = s l or (Ins) 1 
[3, Q> etc. We show that for such wave functions the contribution of B vanishes. The 
contribution of integral B also vanishes for the exact wave function. Thus, in calculations 
of the integral J f(s,t,u)dr when the space considered is infinite, one may ignore the 
contribution of the component B by choosing an appropriate decay of the approximate 
wave function. Nevertheless, we reiterate that from a rigorous mathematical perspective, 
the integral f f(s,t,u)dr is composed of two components A and B. 



APPENDIX 

For the wave function of Eq.(21), the contribution to the total energy from integral A is 

E A = 2n 2 f R ds f'dt I" du[u(s 2 -t 2 )(^) 2 - 4ZsuC 2 e- 2as + (s 2 ~t 2 )C 2 e- 2as )], (A.l) 
Jo Jo Jt as 

where the first, second, and third terms correspond respectively to the kinetic, external, and 
electron-interaction energies. (The surface contribution to the kinetic energy is discussed 
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below). Thus, 

tt 2 C 2 1 
E A = — — [-(27 -8a) - — {e- 2aR [-405 - 30(-4 + 27i?)a-30 J R(-8 + 27i?)a 2 
8a 5 15 

-60,R 2 (-4 + 9R)a 3 + 10(16 - 27R)R 3 a A + 80i?V + 32R 5 a 6 ]}]. (A.2) 

The contribution from integral B is 

»2R rfli-s r s 



f 2R f 2R-s f s dj 

E B = 2tt 2 / ds dt du[u{s 2 -t 2 ){^-f ~4ZsuC 2 e 

Jr Jo Jt vs 



2 -las 



+ (s 2 -t 2 )C 2 e- 2as )} (A.3) 

7r 2 r 2 

= -^e~ 4aR [-21 - 4(-2 + 2lR)a + 8(4 - 17R)Ra 2 - 32i? 2 (-2 + 3R)a 3 

+Q4R 3 a 4 + 32R 4 a b ] 
7r 2 C 2 

+ Y^F e ~ 2ai? t 315 + 30 (~ 4 + 2LR ) a + 30i?(-8 + 5R)a 2 - 60i? 2 (4 + 9i?)a 3 
+10(16 - 27 R)R 3 a 4 + 80R 4 a 5 + 32R 5 a 6 }. (A.4) 

The kinetic energy is 

T = - J iN\i)dr = J (ViVO • (Vi^)dr - ^ V • (^Vi^)dr. (A.5) 

The first term of Eq.(A.5) corresponds to the first terms in Eqs.(A.l) and (A.3) for the total 
energy for integrals A and B, respectively. This contribution is 

T A + T B = a 2 . (A.6) 

The surface contribution to the kinetic energy is 

T s = - y Vi • (V-Vi^dr = " / Vi • (e' ari V ie - ari )rf ri 

= - y (e~ ari Vie- ari ) • dSi (A.7) 

4r> 4 /? 2 p" 2ajR 

In spherical polar coordinates, the kinetic energy T in a volume of radius R is 

-R r R 



so that 



T = - y ^V 2 ^r = 167r 2 C 2 y y e - 2a(ri+r2) (a 2 - ^)r 2 ir 2 2 dndr 2 , (A.9) 

r _ ^ 2 [l+e- 2 ^(-l-2^ + 2a 2 ^ 2 )] 

[1 - e - 2 ^(l + 2a J R + 2a 2 J R 2 )] ' { } 



The sum of Eq.s(A.6) and (A.8) equals Eq.(A.lO). 
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